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Recently two methods have been developed for the quan- 
tization of the electromagnetic field in general dispersing and 
absorbing linear dielectrics. The first is based upon the intro- 
duction of a quantum Langevin current in Maxwell's equa- 
tions [T. Gruner and D.-G. Welsch, Phys. Rev. A 53, 
1818 (1996); Ho Trung Dung, L. Knoll, and D.-G. Welsch, 
Phys. Rev. A 57, 3931 (1998); S. Scheel, L. Knoll, and D.- 
G. Welsch, Phys. Rev. A 58, 700 (1998)], whereas the second 
makes use of a set of auxiliary fields, followed by a canoni- 
cal quantization procedure [A. Tip, Phys. Rev. A 57, 4818 
(1998)]. We show that both approaches are equivalent. 



I. INTRODUCTION 

With the advent of modern optical materials, such as 
optical fibers and photonic crystals, the problem of quan- 
tization of the electromagnetic field in dielectrics has be- 
come an important subject and much activity has been 
taking place in this field. Quantization is required to de- 
scribe the decay of embedded atoms (for specific cases, 
see the Casimir effect and other nonclassical 

phenomena such as the propagation of entangled states 
through dielectrics We also mention the generation 
of X-ray transition radiation by fast electrons traveling 
through layered dielectrics Q . 

In the case of linear conservative dielectrics, quanti- 
zation is well-known for systems, where the permittivity 
(electric permeability) e is a real constant Q or a real 
function of space, e = e(r) [Hjoj. Nonhnear dielectrics 
are discussed in pO| ] . For dispersing and truly absorbing 
media the situation is more complicated, because the per- 
mittivity is a complex function of frequency and varies 
with space in general, i.e., e = £(r,a;). Progress in this 
field has been fairly recent. 

Two basic approaches can be distinguished. The first is 
based on the Hopfield model of a bulk dielectric . The 
quantized electromagnetic field is coupled to a material 
system described by a harmonic oscillator model, and the 
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Hamiltonian of the total system is diagonalized ||T^. A 
drawback is that it becomes rather cumbersome if spatial 
inhomogeneities are present jlSj. Also the identification 
of the permittivity is not trivial [ p^jp^ . 

The second approach starts off from the classical phe- 
nomenological Maxwell equations, featuring a general 
spatially inhomogeneous, complex, frequency-dependent 
permittivity e{T,Lo) satisfying the Kramers-Kronig rela- 
tions. It has the advantage that the really measured val- 
ues of the permittivity can be used for the theoretical de- 
scription of quantized light in media. For example, in the 
case of photonic crystals made up from dielectric objects 
(scatterers) in a conservative, homogeneous background 
(such as vacuum), e(r, u) is known at the outset. Absorp- 
tion is often undesirable, in particular if one is interested 
in band-gap phenomena. In fact absorption prohibits the 
formation of the latter |jl^ . Band-gap photonic crystals 
offer many interesting technological applications (l7| but 
require a large dielectric contrast between the scatter- 
ers and background. This can be accomplished by us- 
ing small metal spheres showing a Drude-type behavior, 
where the real part of the permittivity can acquire large 
negative values |jl^. However, such systems are always 
somewhat absorbing. On the positive side, absorption 
may be advantageous in the case of transition radiation, 
where it can be used to suppress undesired frequencies 

There have been two concepts of quantization of the 
phenomenological Maxwell field for general dispersing 
and absorbing linear dielectrics. The first (referred 
to as LN concept) is based upon the introduction of 
Langevin noise current (and charge) densities, as dictated 
by the fluctuation-dissipation theorem, into the classi- 
cal Maxwell equations, which can then be transferred 
to quantum theory by conversion of the electromagnetic 
field quantities into operators. After some earlier work 
]l^ , p^ , restricted to specific simple geometries, a general 
formalism was put forward by some of us M-Vt- In 
this scheme the dyadic Green's function associated with 
the classical (inhomogeneous) Helmholtz equation plays 
a prominent role. Its properties come into play by deriv- 
ing the equal-time commutators for the fields, given those 
of the noise current operator. In Ref. [|l| the case of a 
planar interface and in Ref. l3,E3] the spontaneous de- 
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cay in a spherical cavity is worked out but more involved 
situations can also be handled. Basically the Green's 
function of the classical problem must be calculated. For 
this, general methods and a variety of specific examples 
are considered in Ref. [|4j. Efficient methods have been 
developed (such as an adaptation of the KKR approach 
of solid state physics) for the photonic crystal case psj . 

The second concept (referred to as AF concept) devel- 
oped by one of us (AT) also starts off from the 
classical phenomenological Maxwell equations. Here, the 
introduction of a set of auxiliary fields (instead of a noise 
current) allows the replacement of Maxwell's equations, 
which feature a time convolution term relating the polar- 
ization to the electric field, by a new set of equations for 
the combined set of electromagnetic and auxiliary fields 
but without time convolutions. For the so extended sys- 
tem, a conserved quantity, bilinear in all fields, general- 
izing the electromagnetic energy, exists. Maxwell's equa- 
tions are retrieved by setting the initial auxiliary fields 
equal to zero. The system can then be quantized and the 
conserved quantity becomes the Hamiltonian. But now 
the initial auxiliary fields, being operators, can no longer 
be set equal to zero. 

Although the two concepts of quantization of the phe- 
nomenological Maxwell field look quite different at first 
glance, they are fully equivalent, as we show in the 
present paper. For this purpose, we first review in Section 
y the basic formulas of the LN concept. After reviewing 



the AF concept, we then prove in Section [II that the AF 



concept precisely leads to the noise current operator in 
the LN concept. Some concluding remarks are given in 



Section IV 



II. THE LANGEVIN NOISE METHOD 

Starting point is the set of the classical macroscopic 
Maxwell equations for the electromagnetic field in an ab- 
sorbing linear dielectric without free charges and currents 



9tD(r,t) =9r X H(r,i), 
dtB{r,t) = -dr X E{r,t), 
ar-D(r,to) =0, 

D(r,t) =£oE(r,i)-KP(r,i), 



(2.1) 
(2.2) 
(2.3) 
(2.4) 
(2.5) 



P(r,t)=£o/ dsx(r,<-s)E(r,s)-HP„(r,i), (2.6) 

Jto 



B(r,i) =/ioH(r,i), 



(2.7) 



where the initial time Iq may be set to to — —oo. Intro- 
ducing the Fourier transform of the electric-field strength 
according to 

/ + 00 
(ia;exp[-iwt]E(r,w) 
-oo 
r+co 

= / dio exp[-iu)t]'E{r, uj) + C.C., (2.8) 
Jo 



and t he F ouri er t ransforms of the other fields accordingly, 
Eqs. (|t|) - (H^ lead to 



dr X H(r,cj) — — icjD(r, (jj), 
dr X E(r,cj) = iuj'B{r,Lu), 
dr-'D{r,Lu)=0, 
dr-Bir,uj)=0, 

D(r, uj) = £oe(r, '^)E(r, lj) + P„(r, lj), 
B(r,w) = ^oH(r,w), 



where 



e{r,uj) l + x{r,uj), 



dt exp[ja;t]x(r, t). 



Note that for absorbing media 

e{Y,uj) = £_R(r,tj) +i£/(r,w), £/(r,w) > 0. 



(2.9) 
(2.10) 
(2.11) 
(2.12) 
(2.13) 
(2.14) 

(2.15) 
(2.16) 



(2.17) 



In the LN concept (2^-|2|j2|, Eqs. (^ - ( pl4| ) [or 
Eqs. ( |2.l| ) ~ (2.7)] are considered as a set of equations 
for the electromagnetic field supplemented with a noise 
polarization Pn(r,a;) Its introduction arises from 

the necessity to fulfill the fluctuation-dissipation theo- 
rem, because macroscopic electrodynamics is a statistical 
theory. In a classical theory the noise term can only be 
dropped in the zero-temperature limit, T — > 0, whereas 
in quantum theory it is always present due to vacuum 
noise. From these arguments, the operator-valued fields 
(indicated with hats) in quantum electrodynamics can be 
regarded as obeying Eqs. ( |2.9| ) - (2.14), 



dr X B(r,u;) 



r£(r,cj)E(r,w)+/xojn(r,c^), (2.18) 



dr X E(r,a;) — iuj'B{T,uj), 

dr ■ £o£(r, cj)E(r, w) = /5„(r, u), 

9r-B(r,c^) = 0, 

D(r, w) = £o£(r, w)E(r, cj) + P„(r, w). 



(2.19) 
(2.20) 
(2.21) 
(2.22) 



where /5n(r, i^j) and jn(r, oj) are the noise charge and cur- 
rent densities, 



/5n(r,Cj) = ~dr ■ Pn(r,cj), 
jn(r, w) = -iwPn(r,u;). 



(2.23) 
(2.24) 



Quantization is accomplished by relating the current to 



bosonic vector fields according to (eo/io - 



jn(l-,w) 



^oc^ V 7r£o 



-£/(r,w)f(r,w) 



= w-\/-^£7(r,cj)f(r,w), 



(2.25) 



f(r,w),ft(r',t^') =5{v-v')5(uj^J)\J, (2.26) 
f(r,c.),f(r',u;')l = fft(r, c.), ft(r', l.')1 =0, (2.27) 
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where U is the unit 3x3 matrix. The fields f(r, w) repre- 
sent the fundamental variables of the overall system. In 
the Heisenberg picture they evolve as 

f (r, oj, t) = exp[-icj(t - t')] f (r, uj, t'), (2.28) 

which is governed by the Hamiltonian 

H = jdvj^ dLohujP[v,u)i{v,uj). (2.29) 

The commutation relations ( 2.26| ) and ( ^.27 ) imply that 



jn(r,u;)jj;(r',u;') 

2 



£/(r, uo)5[v ~ r')5{uj - uj')\} 



= eoe/(r, uj)S{r — r )S{iLj ~ lu )U. 



(2.30) 



From Eqs. ( ^.18D and ( pJ9| ) it follows that E(r, uj) sat- 
isfies the equation 

9r X 9r X E(r, w) ^£(r, a;)E(r, w) 

c/ 

^2 

= [Ho - ^£(r, cj)U] • E(r, w) = ic^Mojn(r, c^) (2.31) 



[Ho = dr X dr X = drdr — d^U) . luvcrsiou of Eq. ( |2.3l| ) 
and the use of Eq. ( 2.19D yields 



E(r, w) = icj^o y dsG(r, s,w) • j„(s,w), 
B(r,w) = (icj)"i(9r X E(r,cj) 

^fJ-odrX / dsG(r,s,w) • jn(s,w). 



(2.32) 



(2.33) 



Here, G is the classical Green function (actually a second- 
rank symmetric tensor) that satisfying the equation 



drdr — 



dt 



-e{v,uj) 



U \ ■ G(r,s,t^) = (5(r-s)U. 

(2.34) 



Note that G corresponds to the operator [drdr — {dr + 
uj^ e{Y , uj) / c^} U]^^, which exists as a bounded Hilbert- 



space operator if 



lim [£{y,uj) - 1] i 0+, 

|r|— >cx3 



(2.35) 



automatically fixing the boundary conditions of G at in- 
finity (cf. Ref. @). 

The electric-field strength operator in the Schrodinger 
picture can then represented in the form of [cf. Eq. (2.8)] 



E(r) = / da;E(r,w) +H.C., 
Jo 



(2.36) 



and the other field operators accordingly. Using the prop- 
erties of G it can be verified that the standard equal-time 
commutation relations of quantum electrodynamics fields 
are fulfilled [^- ^ , ^ . Since the latter do not depend 
on e{v,uj), the case that £i{v^uj) (approximately) van- 
ishes in a certain w-interval can be handled by means 
of a limiting procedure. It is worth noting that the LN 
method has the advantage that arbitrary inhomogeneous, 
anisotropic, amplifying or magnetic matter can easily be 
included in the formalism lE8|. 



III. THE AUXILIARY FIELD METHOD AND ITS 
RELATION TO THE LANGEVIN NOISE 
METHOD 

A. Classical formalism 

The AF method |27j starts from the zero-temperature 
classical Maxwell equations (?„ — 0) and complements 
them with appropriately chosen auxiliary fields. In order 
to facilitate a comparison with the LN method, we shall 
use a setup where only Fourier components for positive 
argument are used and in addition we shall use a different 
gauge for the fields. We assume that x(r, t = Q) = Q, 
which can be verified from linear response theory. It 
excludes instantaneous surges at the initial time. Then, 
with x'(r,t) =9tx(r,t), 

atE(r, t) = c^dr X B(r, t) ~ f ds x'(r, t - s)E(r, s) 

J —OG 

= c^drxB{r,t)-J{r,t), (3.1) 

where J(r,t) = dtP{r,t) is the polarization current den- 
sity. Since x(i";0) = we have [the factor 2 arises from 
changing the range of the A-integral from R in Ref. |2^] 

to [0,oo)] 



(3.2) 



PCX) 

X'{r,t)^2 dX cos{Xt)i^{r,X), (3.3) 

where i/(r, A) > for absorbing systems considered here. 
Note that 



£/(r,^) = ^J^(r,A) 



(A > 0). Next we define 



Fi(r,t) = V^E(r,i), 
F3(r,t) = _LB(r,t) 



(3.4) 

(3.5) 
(3.6) 



and introduce the auxiliary fields 
F2(r,A,t) =-V^a(r,A) 

X / ds sin A(t — s)E(r, s), 

J — oo 



(3.7) 
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F4(r,A,t) 



eocr(r, A) 



F;,(r,A,t) =exp(~iAi)F;,(r, A) 



/ ds cos A(t — s)E(r, s), 



where 

2i/(r, A) = cr(r, A)^ cr(r,A)>0, 
and note that 

F2(r, A, -oo) F4(r, A, -oo) = 0. 
It can be proved that the set of equations 
atFi(r,t) =c9r X F3(r,t) 



(3.8) 

(3.9) 
(3.10) 



+ 



/•OO 

/ dAcr(r, A)F4(r, A,t), (3.11) 
Jo 



5tF2(r,A,i) = AF4(r,A,t), 

5tF3(r,t) = -carXFi(r,<), 

atF4(r, A, t) = -AF2(r, A, t) - aiy, A)Fi (r, t) 



(3.12) 
(3.13) 
(3.14) 



together with the initial conditions ( 3.1Cl| ) is equivalent 
to Maxwell's equations, and the quantity 

£ = \ldv [E(r,t)2+B(r,<)2] 
+ i/dr/g"dA [F2(r,A,i)2+F4(r,A,i)2] (3.15) 

is conserved in time. Note that £ coincides with the 
electromagnetic energy for vanishing x- 



Our aim is to find a quantized ve rsion of Eqs. (3.11) 
- (3.14). Since the initial condition ( |3.10 ) then loses its 
meaning, we now drop it. Setting 



Fo(r, A, t) = F4(r, A, t) - iF2(r, A, t), (3.16) 
we have 

9tFo(r, A, t) = AFo(r, A, t) - a(r, A)Fi (r, i). (3.17) 
Its solution can be written as 
exp(iAi) Fo(r, A, t) 
= F(,(r,A)-cr(r,A) / ds exp(as) Fi(r, s). (3.18) 

•/ — oo 

Since the second term on the right hand side vanishes as 
t — > — oo, we obtain (for the limit, see the Appendix) 

Fo(r,A)= lim exp(at) Fo(r, A, i). (3.19) 

On the other hand, setting ^ = in Eq. ( ^.18| ) yields 
F;,(r,A) =Fo(r,A,0) 
+ CT(r,A) / ds exp(as)Fi(r,s). (3.20) 

J — OO 

With 



Eqs. (3.11) - (3.14) are then replaced by 

atD(r,i)=arXH(r,i), 
atB(r,t) = -^r X E(r,i), 
5tF^(r,A,t) = AF;(r,A,i), 
5tF^(r,A,i) = -AF'2(r,A,i), 



(3.21) 



(3.22) 
(3.23) 
(3.24) 
(3.25) 



, / dsx(r,t~s)E(r,s)+P'(r,t). 

./ — oo 



where 

D(r,i) =eoE(r,t) 

+ £o 

Here, 

P'(r,i) = y" dAA-V(r,A) 

X [cos(At) F^(r, A) + sin(At) F;(r, A)] 
can be regarded as being the noise polarization, with 

J'(r,i)=a*P'(r,<) = y" dAa(A) 



(3.26) 



(3.27) 



[sin(Ai) F'2(r, A) - cos(Ai) F:i(r, A)] 



(3.28) 



being the assoc iated noise current density. In Eqs. (3.27) 
and Eqs. ( |3.28D we have set Y'Q = Y\-i'¥'^ with F^ and F^ 
real. Note that the equations of motion for the primed 
auxiliary fields are decoupled from those of the electro- 
magnetic fields. 



B. Hamilton formalism and quantization 



We can interpret Eqs. ( 3.22| ) - ( 3.25| ) as a set of equa- 
tions suitable for transferring to quantum theory. The 
equations of motion for the F'-fields describe harmonic 
motions and are readily quantized, thus leading to a 
quantum noise contribution in the field equations. How- 
ever, there actually exists a Hamiltonian formalism, gen- 
erating the full set of field equations, which can then be 
:i uant ized ||2^,^. The basic equations are Eqs. ( 3.11 ) - 
.14|), which can be written in the compact form of 



aF = NF = 







(3.29) 



where Eg consists of Fi and F2 and E^ of F3 and F4. 

In order to show the equivalence to the LN method, we 
adopt a generalization of the temporal or Weyl gauge (the 
T- gauge in Ref. H^) instead of the generahzed Coulomb 
gauge in Ref Wh- Thus we set 



F 

F 

A 7 



F3 
F4 



-dt 



^2 



(3.30) 
(3.31) 
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which imphes that 



(3.32) 

(3.33) 
(3.34) 
(3.35) 



B = 5r X £(, ^/^li 

F4 - (7^1 + A|2- 

-1/2 

I 

grangian generating the equations of motion is 
where the (real) inner product is defined according to 



Obviously, A = Eq is the vector potential. The La- 

of mo 

5(Hele,Ce>, (3-36) 



(f , g) = £0 / dr 



fi(r) • gi(r) 



dAf2(r,A)-g2(r,A) 



The canonical momentum is 



7ri(r) 
7r2(r, A) 



(3.37) 



(3.38) 



thus E = — £q ^'^TTi, and the Hamiltonian becomes 



-1/2 
H- 1 



2 \" ei e/ > 2 



(3.39) 



Here (cf. Ref. 



He = N:„,N™e: 



'c2Ho + x'(r,0) J^dXXa{r,X)- 
A(7(r, A) A^ 



(3.40) 

where Hq = drdr — d^U is a non- negative selfad- 
joint operator in Ti^ = Hi ® 7^2, = ^^(R^, dr; R^), 



, drdA; 



Note that the ofF-diagonal 



terms in Eq. ( 3.4C ) describe the coupling of the electro- 
magnetic and auxiliary fields, but that there is also a 
shift term x'(r, 0), and it can be verified that H equals 
the conserved quantity £. 

Quantization is achieved by setting 

[^i(r),7^i(r')] -iM(r-r')U, (3.41) 
[^2(r, A), 7t2(r', A')] = *M(r - r')S{X - A')U, (3.42) 

all other commutators being zero. The Hcisenbc rg eq ua- 
tions of motion are given according to Eqs. ( |3.1l|) ~ 
( PH ) [or Eq. (|^)]. Defining Fo(r,A,i) according 
to E g. (|3.16 ), it again satisfies the (operat or- va lued) 



Eq. ( p.lTI ) with the solution according to Eq. ( 3.18 ): 
Fo(r,A,t) =exp[-^A^]F^,(r, A) 
-cr(r, A) / dsexp[-iX{t - s)]Fi{r,s), (3.43) 



where, according to Eq. ( ^.ig| ) together with Eq. ( |3.16| ), 

F;,(r,A) =Fl(r,A)-*F'2(r,A), (3.44) 

with F'y ^(r, A, t) and F2 4(r, A) selfadjoint. Insertion 
of Eq. (3.43) in the equ ation s of motion results into 
the opera tor- valued Eqs. ( 3.22| ) - ( 3.27 ). According to 
Eq. (3.28|), the operator of the noise current density reads 



J'(r,i) J dX(T{X) 

X sin( At ) F2 (i", A) - cos( At ) F4 (r , A) 



(3.45) 



Using Eqs. (3.41) and (3.42), we find the (equal-time) 
commutation relations 

[F;(r, A), F^(r', A')] = -inX6{r - r')6{X - A')U (3.46) 

and 



[j'(r,A),J't(r',A')] 

fe\2 

= — £o£/(r, X)6{r - r')S{X - A')U. 



(3.47) 



Thus, we identify J'(r, A) with the Langevin noise current 
m Section |l|. 



J'(r,A)=j„(r,A). 



(3.48) 



With this choice, Eq. ( p.47| ) exactly equals Eq. ( p.30[ ). 
Hence, the AF formalism is equivalent with the LN for- 
malism. Introducing creation and annihilation operators 
according to 

Fi(r, A) = z{hX/2)'/^ [b(r. A) - bt(r. A)] , (3.49) 
F;(r, A) = (/iA/2)i/2 [b(r, A) + bt(r. A)] , (3.50) 

we have 

[b(r,A),bt(r',A')] = S{r - r')5{X - X')U, (3.51) 
and hence 

f (r. A) = -b(r. A) = -{2nX)-'/^F'o{r, A). (3.52) 
IV. DISCUSSION 



In the LN formalism, the basic ingredient is the iden- 
tification of the (usually discarded) noise polarization in 
Eq. (2^) or Eq. ( ^.13 ) as the fundamental field variable 
of the theory from which all properties of th e elec tr omag - 
netic fiel d can be derived by means of Eqs. (^ ), (E2I, 
( ^.32D , ( |2.33| ), and Maxwell's equations (^) - ( p.l4| . 
For any temperature, the fluctuation-dissipation theorem 
is then satisfied, and the classical, statistical noise po- 
larization can be regarded as being an operator-valued 
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quantity in quantum theory. Then one can show that 
the correct (equal-time) QED commutation relations are 
satisfied. 

The AF formalism starts from Maxwell's equations 
without the noise polarization. Instead, auxiliary fields 
are introduced whose equations of motion eventually de- 
couple from those of the electromagnetic fields leaving 
behind a source term in Ampere's law. The formalism 
can then be cast into a Hamiltonian form, which, upon 
quantization, features a noise current with the same com- 
mutation properties as in the LN formalism. For ease 
of comparison, a generalized temporal gauge has been 
adopted, but the actual choice does of course not affect 
the equal-time commutator [E(r), BfrQ]. 



It should be pointed out that Eq. ( p.3lD , which plays 
a crucial role in the LN concept, can be related to the 
eigenvalue problem asso ciated wi th He in the AF formal- 
ism. Indeed, from Eqs. ( |3.1lD - ( |3.14D we have 



(4.1) 



Thus, in the stationary solution Fe(t) — exp[—iujt]F , 
Fe{u!) solves the eigenvalue problem 



where the first of these equations, 

[c^Ho + x'm] Fi(^) 

dAAcr(A)F2(A,w) ^uj^Fi{uj), 



(4.2) 



(4.3) 



corresponds to Eq. (|2.3l| ), given the relation ( p. 191) be- 
tween Fq and Fq. As shown in the Appendix, this rela- 
tion has a precise scattering theoretical background. 
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APPENDIX: THE RELATION BETWEEN Fo 
AND Fo 



The existence of the limit 



F;,(r, A)= lim exp(iAt) Fo(r, A, i) 

t — ^ — OO 



(Al) 



suggests that Fo(r, A,f) is asymptotically free, i.e., its 
motion becomes decoupled from that of the electromag- 
netic fields as t — > — oo. Such an asymptotic behavior 



can be studied more precisely in terms of M0ller wave 
operators, as we shall now briefly discuss. We write 



V 






-a(r,A) 



cdr X dX a(r. A) • • • \ 



A 





No + Ni, 



where 



/ 








C(9rX 


"A 















— cdrX 











V 





-A 





0/ 



(A2) 



(A3) 



Ni = 



Let 



/ J^dXa{r,X)- 





\ -cr(r,A) 



/O 
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V 1 



(A4) 



(A5) 



be the projector upon the auxiliary fields. Then it can 
be shown by standard methods that the M0ller operators 



±= lim exp(-Ni) exp(Noi) Pa 

»iboo 



(A6) 



exist in the strong sense (for the underlying Hilbert space 
and other details, cf. |2^). A quite more elaborate analy- 
sis shows that the adjoints rij_ can also be given as strong 
limits, i.e.. 



= lim Pa 



;exp(-Not) exp(Nt) 



exist, implying that 
P 



F{t) exp(NoO f^iF(O), 



(A7) 



(A8) 



i.e. the motion of the auxiliary fields becomes decoupled 
from that of the electromagnetic ones for large times. 
Physically, this can be understood by observing that for 
a finite dielectric the auxiliary fields do not propagate 
(they are confined to the dielectric), whereas the elec- 
tromagnetic fields propagate away. This gives a ri gorou s 
underpinning of the existence of the limits in Eq. ( p. 19 ). 
We can arrive at Eq. (3.20), starting from Eq. (A8), not- 
ing that 

niFiO) = PauxF(O) 

dt Paux exp(-Not) (N- No) exp(Ni) F(0) 

PauxF(0)-f/ diPauxexp(-NoONiF(i), (A9) 
Jo 



6 



and working things out, which leads to 



ds 



in*_F)2\ ^ /F2(0) 
sin(As) 



cos(As~ 



a(A)Fi(s). 



(AlO) 



Following the procedure gi ven a bove, one can show an 
equivalent expression to Eq. ( 3.20|) valid in the quantum 
case. 
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